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^vq . Numerical computations of transport coefficients at low temperatures are pre- 

sented for shapes typically encountered in nuclear fission. The influence of quan- 
tum effects of the nucleonic degrees of freedom is examined, with pair correlations 
included. Consequences for global collective motion are studied for the case of the 
decay rate. The range of temperatures is specified above which this motion may 
' be described as a quantal diffusion process. 
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In the past decade much progress has been made in the understanding of nuclear 
transport phenomena in the regime of not too low temperatures, say between 1 and 
5 MeV (with ks = 1). Such a situation is reached if two heavy ions collide at an 
energy above the Coulomb barrier, but where the excess energy per particle still is small 
compared to the Fermi energy. In this regime the dynamics of the composite system 
may be parameterized in terms of shape variables. Of particular interest is the outgoing 
channel which is dominated by fission and the emission of light particles and 7's. It has 
been possible experimentally to deduce valid information on the time scale of collective 
motion [I]], and, hence, on the size of nuclear dissipation. These experiments suggest 
collective motion to be over-damped, possibly providing an answer to the question raised 
by Kramers as early as 1940 in his seminal paper ||, namely whether nuclear friction 
is "abnormally small or abnormally large". Nowadays such processes are described 
theoretically in terms of the Langevin equation [[| , which is understood to be equivalent 
to Kramers' original equation (of Fokker Planck type) for the density in collective phase 
space. 
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On general grounds, it may be anticipated that the magnitude of nuclear dissipation 
will vary with excitation. Indeed, there are experimental indications for such a 
conjecture. At small thermal excitations the dynamics is governed by the (real) mean 
field for which there is no room for damping of slow collective motion like fission. At 
larger T coupling to more complicated configurations sets in, which causes transfer of 
energy from the collective degrees of freedom to the nucleonic ones Xj. Within the 
linear response approach || the effects of this coupling are accounted for by dressing 
the single particle energies by complex self-energies depending both on frequency and T. 
Approximating its imaginary part by a constant proportional to T 2 , friction will again 
decrease with T, once the microscopic damping becomes so large that one may speak 
of "collision dominance". At intermediate temperatures there might be the intricate 
contribution to friction from the "heat pole", which has been seen to be large for non- 
ergodic systems ||, and which has a dramatic influence on the T-dependence. 

In the present Letter we want to focus on very low excitations, say in the range below 
T « 1 MeV. This regime not only is of great practical importance, as for the production 
of super-heavy elements 0, for instance, it is also of great theoretical interest. First 
of all, there is little doubt that in this domain quantum effects dominate nucleonic 
dynamics, and the transport coefficients will strongly be influenced by shell effects and 
pair correlations. One even must expect quantum features to be present for collective 
motion, for instance as corrections to Kramers formula for the decay rate ||. Often 
quantal approaches are based on the functional integral method applied to simplified 
Hamiltonians of the Caldeira-Leggett type (for a review see ||). There, the bath degrees 
of freedom Xi are represented by a set of oscillators of fixed frequencies, with a bilinear 
coupling between the Xi and the collective variable Q. The decay rate is calculated for 
imaginary time propagation. Both features hardly can be taken over to nuclear fission. 
First of all, the simplest constraint to warrant self-consistency between the mean field 
and the shape of the nuclear density requires the former to change with Q. This aspect 
alone makes it very difficult to work with a (pre-fixed) Hamiltonian for the total system 
of all degrees of freedom. Moreover, the temperature, which one may define [ID] for 



the fast degrees of freedom (supposedly be given by the "nucleonic" ones), is subject to 
changes with Q as well as with time. The latter feature occurs because of the evaporation 
of particles mentioned before. 

For these reasons a formulation with real time propagation is much more appropriate. 
This has been achieved by a suitable application of linear response theory on the basis of 
a locally harmonic approximation (LHA) (for a review see||). One exploits the concept 
of propagators which move the system forward in collective phase space by small time 
steps. As the individual ones only cover small areas they may be represented by (multi- 
dimensional) Gaussians. The latter satisfy an equation of motion whose structure is 
similar to that of Kramers, with only the diffusive terms being modified to account for 
quantum effects. 

The following study is based on numerical calculations of transport coefficients for 
average motion, namely friction 7, inertia M and local stiffness C, more precisely of 
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those ratios which determine transport in phase space, 
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Their knowledge will allow us to examine implications on the diffusion coefficients and, 
hence, on transport processes like fission. It would be most desirable that such informa- 
tion be used as input for computational codes to solve for Fokker-Planck or Langevin 
equations. In this way one would be able to examine in more detail the role of shell 
effects, which are known to produce structure not only in the static energy but in the 
inertial and frictional forces as well. To simplify matters, for the present purpose we will 
look at the more schematic case where the system's energy exhibits just one minimum 
and one barrier at Q a and Q b , respectively. The stiffnesses and the barrier height are 
found from a Strutinsky calculation of the free energy. Finer details of shell effects are 
removed both from the potential as well as from the transport coefficients by applying 
a suitable smoothing over a small region of the collective variable around Q a and Q b . 

Suppose we may at first discard any quantum effects in the collective degrees of 
freedom, which amounts to look at the "high temperature limit", for which Kramers' 
equation applies ||. The temperatures we have in mind are always small compared to 
the barrier height, T <C E b . The decay rate Rk then shows the following behavior. 
For given T, but as function of the r\ b (at the barrier), the RKij] b) T) increases first, to 
decrease after it has reached a maximal value (see e.g. |IT]). The decreasing branch is 
represented well by Kramers' "high viscosity limit" 
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Rk = 2^ [y 1 + %2 ~ Vb ) ex ^~ Eb ^ ^ 

which is valid for r\ b > T / (2E b ) (see e.g.||). If blindly extended down to r] b = this form 
Rk reaches a value typical for a simple transition state model, ttst ~ r {Vb — 0) (Bohr- 
Wheeler formula). Rather, for very small rj b , one ought to apply Kramers' "low viscosity 
limit", given by R l K J ^ = T h ^(E b /T) exp(— E b /T). For nuclear physics the latter has not 
played any role yet, as rj is believed to lie above the limit given below (0). According 
to H H, pi this should be the case at temperatures above 1 MeV. Moreover, the w 



does not change much with T. In |12| a value of about 1 MeV/h was found both at the 
potential minimum as well as at the saddle. It so turns out that this feature is more or 
less recovered even at smaller T, say within an accuracy of the order of 20%, which may 
be good enough for the following discussion. 

More drastic modifications are expected, and indeed seen, for dissipation. To study 
this behavior, the 7, T 7 and rj have been calculated on the basis of the same deformed 
shell model as in |E| but with pair correlations included. The transformation from 



independent particles to quasi-particles of BCS-type is standard. For our purpose the 
common procedure does not suffice, however. As mentioned previously, a decent and 
sensible description of nuclear dissipation needs to account for "collisions". At low 
thermal excitations their effects, too, will strongly be influenced by pair correlations. 
Look at the extreme case of zero temperature and take an even-even system for the sake 
of simplicity. Then there will be no quasi-particle states within twice the gap energy 
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2A. Hence, the imaginary part of the self-energy T(huj, A, T = 0) must be zero at least 
within such a range of frequencies hut. Hence, at T = friction will strictly vanish for any 
collective motion whose frequency u — w lies in that range, i.e. Two < 2A. Extrapolating 
from the case of T = to finite T we should expect the function 7 = j(T) to have a 
step like behavior. This dependence then goes through to that of T 7 and rj, albeit the 
inertia, too, is influenced by pairing. 
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Figure 1: The damping factor rj — 7/ (2<J M\C\) at the potential minimum (top) and at 
the fission barrier (bottom) of 224 Th; the dashed curves are for A = 0. The long-dashed 
line in the bottom part shows the function T/2Ej,. 



Let us demonstrate these features on the basis of numerical calculations performed 
for the example of 224 Th. For details we have to refer to [|14[], but we may mention that 



the T(fwj, A, T) has been calculated along the lines suggested in |TJJ. In Fig.l we display 
the 7y(T)'s at the minimum and the barrier. They have been obtained for a A = A(T) 
as determined by the gap equation. Unfortunately, so far it has not been possible to 
calculate the underlying response function xif^) i n full gl° r y- Rather, when evaluating 
the necessary folding integrals over frequency the correct width had to be approximated 
by a constant, calculated at the Fermi energy \i: T(hu, A, T) = T(huj = /x, A, T). Indeed, 
in this regime of small uj and for T < T pair where pair correlations disappear, such an 
estimate may be considered to represent the correct width well enough to allow for a 
general analysis. Evidently, the values of 77's obtained for A / clearly fall well below 
those of the unpaired case, shown here by the dashed lines. 

The most important features exhibited in Fig.l may be summarized as follows, to- 
gether with the consequences for Kramers' decay rate: (1) The step like function men- 
tioned before is clearly visible. (2) Below T « T pa i T ~ 0.5 MeV the effective damping 
rate rj is smaller than about 0.1. (3) As seen in Fig.l, 77 may fall below T/2Eb such 
that formula (^) no longer applies. (4) Up to T ~ 1 MeV rj stays below ~ 0.2 at the 
minimum and below ~ 0.3 at the barrier. The latter value implies that the rate may be 
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approximated fairly well by the transition state value ttst ~ Rk(v = 0), see (fj). (5) 
These values of 77 are much smaller than those one gets within "macroscopic models", 
say in terms of a combination of wall friction with the stiffness and inertia of the liquid 
drop model (with irrotational flow). 
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Figure 2: The ratio of the Kramers correction factor fx = y 1 + Vb~Vb to its macroscopic 
counterpart as function of temperature. The quantity /q/r-/ j™ acro accounts for the 
quantum correction factor /q, see eq.(||). 



In Fig. 2 we plot the ratio of Kramers' pre-factor f K = (Jl + rfc — rjb) obtained from 



our results to that of the macroscopic limit just described flq| . It is seen that the latter 
underestimates the decay rate by about a factor of 10. Quantum corrections will increase 
this deviation further, indicated here by the dashed curve. 

Within the LHA, these quantum corrections come in through the diffusion coeffi- 
cients, as given by the fluctuation dissipation theorem ||. It is only at temperatures 
above 2 MeV that one may safely assume the classic Einstein relation D pp = 7T to be 
valid 0. In the general case, in addition to the D pp there is a cross term D qp , both of 
which depend in non-linear way on combinations of M, 7 and C, or on the parameters 
introduced in (|]). The diffusion coefficients behave very differently for stable and for 
unstable modes. To demonstrate this feature let us look at the limit of small dissipation 
77 <C 1. To lowest order in 7 one gets 



D qp = D pp = 7T* with T*{w) 



hw 



coth 




(3) 



with w = \jCjM = \zu\ for C > and w = i\zu\ for C < 0. The form (^) may be said 
to represent the correct behavior fairly well below r) ~ 0.1 (see Fig.3.4.2 of ||). From 
the results shown above for 77 (T), one may thus argue the relation @ to be acceptable 
for temperatures below T pair , whereas deviations must be expected for T > T pair . For 
C < and weak friction the diffusion coefficient D pp falls below the values given by the 
Einstein relation. It quickly drops to zero at a critical temperature T c , below which the 
D pp would become negative and the diffusion equation would loose its meaning. The 
value of this T c decreases with increasing 77, such that the form given in @ delivers an 
upper limit and we may write 



Tc < T c {r] = 0) = hm b /ir 



T, < T, 



pair 



(4) 
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The statement on the right is reached assuming the h,w b to be of the order of 1 MeV 
and taking the value for T pair , as reported above, together with the fact that below T pair 
the damping rate r\ falls below 0.1. 

Commonly, the quantum corrections to Kramers decay rate are expressed by a factor 
Jq appearing in the correct rate as R = /qRk (see e.g. |§). As shown in ||, this form 
may also be obtained within the LHA. This derivation is based on the assumption 
that in the neighborhood of the potential minimum friction is large enough to ensure 
sufficient relaxation inside the well. The same assumption is behind Kramers' "high 
viscosity limit", upon which we just have convinced ourselves to be given in the range of 
temperature at T pair and above. Moreover, this assumption turns out to be necessary also 
when the problem is formulated and solved with path integrals in real time propagation 
(see 0). 

The /q can be expressed by a ratio of two partition sums: fq = \Z b \/Z a , where the 
one associated to the barrier has to be defined by analytic continuation. According to 
the Z of a damped oscillator can be calculated from the equilibrium fluctuations 



of momentum and coordinate. Hence, within the LHA it might be expressed by the dif- 
fusion coefficients. Unfortunately, for 7 ^ a calculation of the momentum fluctuation 
requires regularization, for instance by introducing a frequency dependent friction coeffi- 
cient (Drude regularization). To get a fairly simple estimate of /q and its T-dependence 
we used the following formula (with hu n = n2irT) 

n=l U n + ^n 1 7 ~ W b 1 

It may be noted in passing, (a) that without (Drude) regularization this formula would 
diverge for r° 7^ and (b) that problems of this type are absent for the Caldeira-Leggett 
approach where the transport coefficients do not change with the collective variable; 
generalizations are possible, though, for instance by introducing variable coefficients 
phenomenologically, see e.g. jBj|. The result of a numerical evaluation of ([5]) within our 



theory is shown in Fig. 2 by the dashed curve. This graph demonstrates several features, 
valid in this range of temperatures: (i) The quantum effects in the collective motion 
may change the decay rate by about 30 % or less, (ii) Already at T = 1 MeV they 
only amount to about 10 %. (iii) More important are the quantum effects of nucleonic 
motion, which are responsible for the deviation of the transport coefficients from the 
macroscopic models. 

Unfortunately, it is not possible to carry the analysis further down to smaller tem- 
peratures. Below T c [recall (|J)] we have seen the LHA break down. It may be said 
that this problem generally appears in real time formulations. Within the functional 
integral approach this has been demonstrated in [17| and traced back to the harmonic 
approximation to the barrier. The problem at stake here is a very severe one for any 
applications of Langevin or Fokker Planck equations to nuclear physics. Both methods 
allow one to account for various non-linear effects, as manifested by variable transport 
coefficients, for instance, but both of them rely on "real time propagation". Moreover, 
practical computer programs exploit locally harmonic approximations, in one way or 
other, such that it is not possible to even define meaningful diffusion coefficients below 
T 

J- c- 
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In this context we should like to mention the possibility of calculating and exploiting 
a Feynman- Vernon functional for global motion on the basis of Random Matrix Theory 
^0| . In principle, this might allow one to study quantum effects, but it is somewhat 



questionable whether this procedure will be applicable to nuclear physics, at those tem- 
peratures where these quantum effects become important which we just mentioned. This 
concern has essentially two reasons — discarding for the moment the very fact that this 
model, too, has difficulties with self-consistency (i) Generally an application of RMT 
ceases to be valid at low excitations, (ii) So far practical applications have been possible 
only to leading order in an expansion in 1/T, actually to that regime of T where friction 
decreases with T. 

Summarizing our results we hope to have been able to exhibit for low energy nuclear 
physics an exciting problem of quantum transport which still is lacking a general solution. 
With respect to the application in nuclear physics itself, it appears to be very difficult 
to describe theoretically processes like the "cold" production of super heavy elements 
without a decent understanding of transport at small temperatures and weak dissipation. 
Whereas low thermal excitations are dictated by experimental conditions, the fact of 
small friction then is a consequence of the quantal nature of nucleonic dynamics in 
a mean field, in particular when pair correlations become important. Of course, to 
obtain more quantitative results, further studies on the microscopic level are needed. 
For instance, it is necessary to understand better the mechanism of "collisions" under 
the presence of pair correlations. Their role on the T-dependence of transport, as well 
as that of the "heat pole" in the larger range of excitations, require further clarification. 
Likewise, it should be very interesting to allow for fluctuations in the gap parameter and 
to examine in which way they might modify fission dynamics. 

The authors like to acknowledge fruitful discussions with J. Ankerhold and N.V. 
Antonenko. This work was supported by the Deutsche Forschungsgemeinschaft. 
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